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Abstract: Given the growing use of Sersic’s (1963, 1968) model for describing the 

stellar distributions in galaxies, and the lack of any single reference that provides the 
various associated mathematical expressions, we have endeavoured to compile such a 
resource here. We present the standard intensity profile, and its various guises such as the 
luminosity, surface—brightness, and aperture—magnitude profile. Expressions to transform 
the effective surface brightness into the mean effective and central surface brightness are 
also given, as is the expression to transform between effective radii and exponential scale- 
lengths. We additionally provide expressions for deriving the ‘concentration’ of an R}^ n 
profile, and two useful equations for the logarithmic slope of the light—profile are given. 
Petrosian radii and fluxes are also derived for a range of Sersic profiles and compared 
with the effective radii and total flux. Similarly, expressions to obtain Kron radii and 
fluxes are presented as a function of the Sersic index n and the number of effective radii 
sampled. Illustrative figures are provided throughout. Finally, the core—Sersic model, 
consisting of an inner power—law and an outer—Sersic function, is presented. 

Keywords: galaxies: structure, galaxies: fundamental parameters, methods: analytical, methods: 
data analysis 


1 Introduction 

Working with the 30-inch Reynolds telescope 1 at what 
was then Australia’s Commonwealth Observatory, and 
today known as Mount Stromlo Observatory, Gerard 
de Vaucouleurs published in 1956 the most extensive 
southern galaxy Atlas of the day. In the following year, 
Jose Luis Sersic commenced work at the 1.54-m tele¬ 
scope at the Astrophysical Station at Bosque Alegre in 
Argentina. His studies from 1957-1966 culminated in 
his 1968 southern-hemisphere galaxy Atlas ‘Galaxias 
Australes’. It too has proven an invaluable contribu¬ 
tion to our understanding of galaxies, evidenced by its 
status as a top 1000 cited astronomy publication. 

In the Introduction of Sersic’s Atlas, it not only 
states the merits for a visual representation of galax¬ 
ies, but, like de Vaucouleurs’, it stresses the necessity 
to go beyond this and obtain quantitative measures of 
the light distribution. This was not mere rhetoric as 
his Atlas consists of two parts, one pictorial in nature 
and the latter quantitative. It is apparent that his gen¬ 
eralisation of de Vaucouleurs’ (1948, 1959) R 1 ^ 4 model 
to an R}^ n model was not merely something he men¬ 
tioned in passing, but something which he felt should 
be done. Indeed, Sersic fitted the R 1 ^" model to every 
(sufficiently large) galaxy in his Atlas. He derived ex¬ 
pressions to compute total (extrapolated) galaxy mag¬ 
nitudes, provided tables of assorted structural parame- 

1 The Reynolds telescope was sadly destroyed in the 2003 

Canberra bush fires. 


ters associated with the R}^ n model, and showed how 
they correlated with galaxy morphological type (his 
Figure 3) and galaxy concentration (his Figure 4, page 
145). Sersic (1963) even provides a prescription to cor¬ 
rect the I? 1 / 71 model parameters for Gaussian seeing 
due to atmospheric and instrumental dispersion. 

It is, however, of interest to note that Sersic’s con¬ 
viction lay in the observation that different galaxies 
possessed differing degrees of an R 1 ^ 4 bulge and an 
R l, ' ] disk component. This mixture of bulge and disk 
components produces a combined surface brightness 
profile with an intermediary form, hence the R}' n model. 

Today, usually when the required resolution is lack¬ 
ing to properly decompose an image into its sepa¬ 
rate bulge and disk components, galaxies are mod¬ 
elled with a single R}^ n profile, just as Sersic proposed 
(e.g., Blanton et al. 2003). While such an approach 
certainly has its merits, we now know that dynami¬ 
cally hot stellar systems themselves posses a range of 
profile shapes that are well described with the R}^ n 
model (e.g., Graham & Guzman 2003, and references 
therein). Detailed studies of well resolved lenticular 
and disk galaxies are routinely fitted with the com¬ 
bination of an exponential-disk plus an R 1 ^" bulge 
(e.g., Andredakis, Peletier, & Balcells 1995; Seigar 
& James 1998; Iodice, D’Onofrio, & Capaccioli 1997, 

1999; Khosroshahi, Wadadekar, & Kembhavi 2000; D’Onofrio 
2001; Graham 2001a; Mollenhoff & Heidt 2001). In 
either case, since the work of Capaccioli in the late 
1980s and in particular Caon, Capaccioli, & D’Onofrio 
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(1993) and D’Onofrio, Capaccioli, & Caon (1994), the 
past decade has seen an explosion in the application of 
the R}^ n model (e.g., Ccllone, Forte, & Geisler 1994; 
Vennik & Richter 1994; Young & Currie 1994, 1995; 
Graham et al. 1996; Karachentseva 1996, Vennik et al. 
1996, to mention just a few early papers), yet no single 
resource exists for the expressions and quantities per¬ 
taining to the R 1 ' n model. Moreover, no one reference 
provides more than a few of the relevant equations, 
and many textbooks still only refer to the R 1 ^ model. 

This (largely review) article intends to provide a 
compendium of equations, numbers, and figures for 
ease of reference. The derivation of these also pro¬ 
vide useful exercises for students. Where appropriate, 
we have endeavoured to cite the first, or at least a 
useful early, reference to any given equation. To the 
best of our knowledge, Figures @ through 1101 . de¬ 
scribing Petrosian indices and Kron magnitudes, have 
never been seen before. A brief reference to where 
readers can find deprojected expressions, and how to 
deal with practical issues such as seeing, is given at the 
end. No attempt has been made here to show the nu¬ 
merous scientific advances engendered via application 
of the R}t n model. 


2 Sersic related quantities 

2.1 The Sersic profile 

Sersic’s (1963; 1968) R 1 ^ n model is most commonly 
expressed as an intensity profile, such that 


I(R) = I e exp 



(1) 


where I e is the intensity at the effective radius R e that 
encloses half of the total light from the model (Caon 
et al. 1993; see also Ciotti 1991, his Equation 1). The 
constant b n is defined in terms of the third and final 
parameter n which describes the ‘shape’ of the light- 
profile, and is given below . 2 

One can integrate Equation CD over a projected 
area A = ivR 2 to obtain the luminosity, L, interior to 
any radius R. This is simply a matter of solving the 
integral 3 


L(< R) = f I(R')2-KR'dR', 
Jo 


which yields, after substituting in x = b n {R/Re) 1 ^", 


where y( 2 n,x) is the incomplete gamma function de¬ 
fined by 


7(2 n,x) = / e t t 2n 1 dt. (3) 


Jo 

Replacing 7(2 n,x) with T(2n) in Equation 0 gives 
one the value of Ltot (Ciotti 1991). 

Thus, the value of b n which we saw in Equation O 
is such that 


r(2n) = 2 7 (2n, b n ), (4) 

where T is the (complete) gamma function (Ciotti 1991) 
Common values of b n are 64 = 7.669 and 61 = 1.678. 
In passing we note a useful property of the V function, 
which is, r(2n) = (2n — 1 )!. 

Analytical expressions which approximate the value 
of b n have been developed. Capaccioli (1989) pro¬ 
vided one of the Erst such approximations such that 
b n = 1.9992n — 0.3271, for 0.5 < n < 10 (see also 
Prugniel & Simien 1997, their equation A3a). Ciotti 
& Bertin (1999) showed b n —* 2n — 1/3 for large val¬ 
ues of n, and in practice this provides a better Et for 
values of n greater than about 8 (see Graham 2001a, 
his Figure 2). Ciotti & Bertin (1999) also provided an 
asymptotic expansion which, for values of n > 0.36, is 
accurate to better than 10~ 4 and the approximation of 
choice. For smaller values of n, MacArthur, Courteau, 
& Holtzman (2003) provide a fourth order polynomial 
which is accurate to better than two parts in 10 3 (see 
their Figure 23). However, the exact value of b n in 
Equation 0 can be solved numerically, and fast codes 
exist to do this, 

For an exponential (n = 1) profile, 99.1% of the 
flux resides within the inner 4 R e (90.8% within the 
inner 4 scale-lengths) and 99.8% of the flux resides 
within the inner 5 R e (96.0% within the inner 5 scale- 
lengths). For an n = 4 profile, 84.7% of the flux resides 
with the inner 4 R e and 88.4%) within the inner 5 R e . 

Multiplying the negative logarithm of the luminos¬ 
ity profile (Equation^ by 2.5 gives the enclosed mag¬ 
nitude profile, known as the “curve of growth”, 


m(< R) = fie 5 log R e — 2.5 log 


2nn 


{b n 


:7(2 n,x) 


(5) 

which asymptotes to the total apparent magnitude mtot 
as R tends to infinity and, consequently, y(2n, a;) —> 
F(2n) (Figured. 

Multiplying the negative logarithm of Equation d 
by 2.5 yields the surface brightness profile (Figured; 
as used in Caon et al. (1993), 


m = Me + [{R/Re) 1/n - 1 ] • (6) 


L{< R) = I e R 2 2nn ° y(2n,a:), (2) 

[On) 

2 It is common for researchers studying dwarf galaxies 
to replace the exponent 1/n with n. In this case, de Vau- 
couleurs’ model would have n = 0.25, rather than 4. 

3 Obviously if one is using a major— or minor-axis pro¬ 

file, rather than the geometric mean ( R = \/ab) profile, an 
ellipticity term will be required. This is trivial to add and 
for the sake of simplicity won’t be included here. The is¬ 
sue of ellipticity gradients is more difficult, and interested 
readers should look at Ferrari et al. (2004). 


2.1.1 Asymptotic behavior for large n 

For large values of n , the Sersic model tends to a 
power-law with slope equal to 5. Substituting e £ = 
z = R/R e into Equation 0, one has 

I(z) ~ exp {-bn. [e t/n - l] } . 

Now, for large n, e 4 /” is small, and so one can use 
e 4 / 71 ~ 1 + t/n. One can also use b n ~ 2n to give 

ln[/(,z)] ~ —bn [t/n] ~ —2 1 ~ —2 In (z). 
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Figure 1: Top panel: Sersic surface brightness 
profiles (Equation 0 for n=0.5, 1, 2, 4, and 10. 
The profiles have been normalised at p e = 20 mag 
arcsec -2 . Bottom panel: Sersic aperture magni¬ 
tude profiles (Equation 0, normalised such that 
the total magnitude equals zero. The dotted line 
is offset by 0.75 mag (a factor of 2 in flux) from 
the total magnitude. 


Thus p(z) = — 2.5 log[/(z)] ~ 51og(z). 

For simplicity, the subscript ‘n’ will be dropped 
from the term b n in what follows. 

2.2 Surface brightness, radial scale, 
and absolute magnitude 

From the value of p e , the ‘effective surface brightness’ 
at R e , and knowing the value of n, one can compute 
both the central surface brightness po and the aver¬ 
age/mean surface brightness (p) e within the effective 
radius. 

At the centre of the profile one has, from Equa¬ 
tion ©, 

p(R = 0) = p 0 = Me - 2.56/ ln(10), (7) 

po = p e - 1.822, n — 1, 
po = Pe — 8.327, n = 4. 


The difference between p e and po is shown in Fig¬ 
ure 0 as a function of the Sersic index n. 

The ‘mean effective surface brightness’, often sim¬ 
ply referred to as the ‘mean surface brightness’, is com¬ 
puted as follows. The average intensity, (7) e , within 
the effective radius is obtained by integrating the in¬ 
tensity over the area A = 'kR/ such that 


(I)e 


JjdA 

A 


I e e b C e~ b(R/Re)1/n 2KRdR 
nRi 



n 

Figure 2: Top panel: Difference between the cen¬ 
tral surface brightness po and the effective surface 
brightness p e as a function of profile shape n. Bot¬ 
tom panel: Difference between p e and the mean 
effective surface brightness {p) e as a function of n. 


Letting x = b(R/R e ) 1//n , one has 
</)e = Je/(n), 

where 


f{n) = 


2 ne° 

6 2 " 


— x 2n— 1 i 

e x ax. 


Now as 6 was chosen such that R e is the radius con¬ 
taining half of the total light, one has 

b r°° b 

/(") = gr J o e-V"- 1 d^^F(2n). (8) 

Thus, 


(p) e = p e - 2.51og[/(n)], (9) 

{p)e = Pe - 0.699, n = 1, 

(p) e = p e — 1.393, n = 4. 

The difference between p e and (p) e is shown in Fig¬ 
ure j2j as a function of the Sersic index n (Caon et al. 
1994; Graham & Colless 1997). 

Substituting equation [15] into Equation [S] one has, 
at R = R e , 

m{< Re) = (p) e — 2.5 log(7rT?e), (10) 

and thus 4 


mtot = (p) e — 2.5 log(27ri?e), (11) 

4 Using empirical measurements within some suitably 
large aperture, one has from simple geometry that (/x) 1/2 = 
^xtot,ap + 2.5 log(27rilJy 2 ). Expressions to correct these ap¬ 
proximate values — due to the missed flux outside of one’s 
chosen aperture — are given in Graham et al. (2005). 
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This expression can be rewritten in terms of the ab¬ 
solute magnitude, Mtot, the effective radius in kpc, 
7? e ,kpc, and the absolute effective surface brightness, 
(fj,) e , aba, (he., the mean effective surface brightness if 
the galaxy was at a distance of lOpc): 


Mtot 


Mtot — (/r) e ,abs 



- 2.51og(27i\R 2 kpc ) 

360 x 60 x 60 \ 2 " 
2 tt x 0.01 ) ’ 


(/r)e,abs - 2.5 log(27r/?e,kpc) - 36.57 


( 12 ) 


The apparent and absolute mean effective surface bright¬ 
nesses are related by the cosmological corrections: 


(M)e.abs = (/r)e - 101og(l + z) - E(z) - K (z), (13) 


2.3 Concentration 

For many years astronomers have had an interest in 
how centrally concentrated a galaxy’s stellar distribu¬ 
tion is (e.g., Morgan 1958, 1959, 1962; Fraser 1972; 
de Vaucouleurs 1977). Sersic (1968) used de Vau- 
couleurs (1956) somewhat subjective size ratio between 
the main region of the galaxy and the apparent maxi¬ 
mum dimension of the galaxy as a measure of concen¬ 
tration. 

Trujillo, Graham & Caon (2001c) defined a useful, 
objective expression for concentration, such that, in 
pixelated form 

Cfl e (a) = JlJ . ( 20 ) 

2^i,j£E(R e ) C 


where z, E(z), and K(z) are the redshift, evolutionary 
correction, and K-correction respectively (e.g., Driver 
et al. 1994 and references therein). 

Another transformation arises from the use of scale- 
lengths h rather than effective radii R e . When the 
R}/ n model is written as 

I{R) = I 0 exp [-{R/h) 1/n ] (14) 

(e.g., Ellis & Perry 1979, their page 362; Davies et al. 
1988), where Iq = I(R = 0), one has 

Io = he\ (15) 

Re = b n h, (16) 

R e = 1.678/i, n = 1, 

Re = 3466/i, n = 4. 

It’s straightforward to show that 

9 5 / fi>\ 1 / n 

• < 17 > 

[i(R) = /ro + 1.086(i?//i), n m 1, 

and 

Ltot — 7r/ 0 /i 2 r(2n + 1). (18) 

Given the small scale-lengths associated with the n = 
4 model, and the practical uncertainties in deriving 
a galaxy’s central brightness, one can appreciate why 
Equation 0 is preferred over Equation fra . 

If one is modelling a two-component spiral galaxy, 
consisting of an exponential disk and an R}M bulge, 
then the bulge-to-disk luminosity ratio is given by the 
expression 


Here, E(R e ) means the isophote which encloses half of 
the total light, and E(aR e ) is the isophote at a radius 
a (0 < a < 1) times R e . This is a flux ratio. For a 
Sersic profile which extends to infinity, 


CrM) 


7(2 n,ba 1/n ) 
7(2 n,b) 


( 21 ) 


This expression is a monotonically increasing function 
of n, and for a = 1/3 it’s values are shown in Fig¬ 
ure 0 . An often unrealized point is that if every el¬ 
liptical galaxy had an R 1 ^ 4 profile then they would all 
have exactly the same degree of concentration. Obser¬ 
vational errors aside, it is only because elliptical galaxy 
light-profiles deviate from the R 1 ^ 4 model that a range 
of concentrations exist. This is true for all objective 
concentration indices in use today. 



B 

D 


zT(2n)e b 

b 2n 


#)(!)' 


(19) 


where h and 7o are respectively the scale-length and 
central intensity of the disk, and R e ,I e , and n describe 
the Sersic bulge profile. Noting that 2nF(2n) = r(2n+ 
1 ) = ( 2 n)!, the above equation can be simplified for 
integer values of 2n. For those who are curious, the 
first term on the right hand side of the equality can be 
seen plotted as a function of n in Graham (2001a). 


Figure 3: The central concentration C/j e (l/3), as 
defined by Trujillo et al. (2001c), is a monotoni¬ 
cally increasing function of the Sersic index n. 


It should be noted that astronomers don’t actu¬ 
ally know where the edges of elliptical galaxies oc¬ 
cur; their light-profiles appear to peeter-out into the 
background noise of the sky. Due to the presence of 
faint, high-redshift galaxies and scattered light, it is 










www.publish.csiro.au/journals/pasa 


5 


not possible to determine the sky-background to an 
infinite degree of accuracy. From an analysis of such 
sky-background noise sources, Dalcanton & Bernstein 
(2000, see also Capaccioli & de Vaucouleurs 1983) de¬ 
termined the limiting surface brightness threshold to 
be p ~ 29.5 B-mag arcsec -2 and p ~ 29 R-mag 
arcsec -2 . Such depths are practically never achieved 
and shallow exposures often fail to include a signifi¬ 
cant fraction of an elliptical galaxy’s light. One would 
therefore like to know how the concentration index 
may vary when different galaxy radial extents are sam¬ 
pled but no effort is made to account for the missed 
galaxy flux. The resultant impact on Cr b and other 
popular concentration indices is addressed in Graham, 
Trujillo, & Caon (2001). 

It was actually, in part, because of the unstable 
nature of the popular concentration indices that Tru¬ 
jillo et al. ( 2001 c) introduced the notably more stable 
index given in Equations m and 1211) . The other 
reason was because the concentration index (7(a) = 

J2i : jeE( a ) -WEi.jesa) Iij ’ where Edenotes some 

inner radius which is a (0 < a < 1 ) times the outer¬ 
most radius which has been normalized to 1 (Okamura, 
Kodaira, & Watanabe 1984; Doi, Fukugita, & Oka¬ 
mura 1993; Abraham et al. 1994), should tend to 1 
for practically every profile that is sampled to a large 
enough radius. It is only because of measurement er¬ 
rors, undersampling, or the presence of truncated pro¬ 
files such as the exponential disks in spiral galaxies, 
that this index deviates from a value of 1 . 


2.4 Profile slopes 




_i_I__1. --— i - --I_ ■ I 

0 2 4 6 8 


R/R e 


Figure 5: The slope a (Equation 1251) is shown as 
a function of normalised radius R/R e for ?r=0.5, 
1, 2, 4, 7, and 10. 


(the inner few hundred parsec) stellar distribution has 
become a quantity of interest. Defining 5 


7 (R) 


-d[log I(R)\ 
dlog R 


( 22 ) 


Rest et al. (2001, their Equation 8 ) used this to mea¬ 
sure the nuclear slopes of ‘core’ and ‘power-law’ galax¬ 
ies. From Equation 0 one can obtain 


7 (R,n) = (b/n)(R/R e ) 1/n . (23) 


This is approximately 2{R/R/) 1 ^ (see section 12 .It . 
Thus, at constant (R/R e ), 7 is a monotonically in¬ 
creasing function of the Sersic index n (Graham et al. 
2003b). 

It turns out Equation is appropriate for the 
so-called ‘power-law’ galaxies which are now known 
to possess Sersic profiles down to their resolution limit 
(Trujillo et al. 2004) and would be better referred to as 
‘Sersic’ galaxies as they do not have power-law profiles. 
A modification is however required for the luminous 
‘core galaxies’, and is described in Section rm 

Another logarithmic slope of interest is that used 
by Gunn & Oke (1975) and Hoessel (1980), and is de¬ 
fined as 


a(R) 


d[lnL(R)] 
d In R 


(24) 


From Equation @1 one has 


Figure 4: The slope of the Sersic profile 7 (Equa¬ 
tion El is shown as a function of profile shape n 
for i?/i? e =0.01, 0.05, 0.1, and 0.2. 


Given HST’s ability to resolve the inner light-profiles 
of nearby galaxies, the slope 7 of a galaxy’s nuclear 


a; d[L(x)] e~ x x 2n 
01 X ’ nL(x) dx ny(2n,x)' 


(25) 


where, as before, x = b(R/R/) 1 ^ (Graham et al. 1996, 
their equation 8 ). 

5 This 7 should not be confused with the incomplete 
gamma function seen in Equation 
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Figures 0 and 0 show how 7 (R) and a(R) vary 
with normalised radius R/R e for a range of different 
profile shapes n. 

2.5 Petrosian index and magnitude 

The Petrosian (1976, his Equation 7) function r](R) is 
given as 


v(R) 


2 tt f* I(R')R'dR' 

7T R 2 I(R) ’ 

L(< R) _ (JU 
7T R 2 I(R) I (R) ' 


(26) 

(27) 


It is the average intensity within some projected ra¬ 
dius R divided by the intensity at that radius. The 
logarithmic expression is written as 


2.5 log [v(R)] = m(-R) - (28) 

and is shown in Figure © for a range of profile shapes 
n. 



R/R e 


Figure 6: The logarithm of the Petrosian func¬ 
tion r](R) (Equation EHl is shown as a function of 
normalised radius R/R e for Sersic profiles having 
n=0.5, 1, 2, 3,... 10. 


where a(R) is given in Equation 1241 Djorgovski & 
Spinrad 1981; Djorgovski, Spinrad& Marr 1984; Sandage 
& Perelmuter 1990, their Section Ha; Kjasrgaard, Jor¬ 
gensen, & Moles 1993). Thus 


r]{x,n) 


2 rry(2n, x ) 
e 


(30) 


The flux within twice the radius Rp when l/g(Rp) = 
0.2 is often used to estimate an object’s flux (e.g., 
Bershady, Jangren, & Conselice 2000; Blanton et al. 
2001), as is the flux within 3Rp when l/r](Rp) = 0.5 
(e.g., Conselice, Gallagher, & Wyse 2002; Conselice 
et al. 2003). How well this works of course depends 
on the shape of the light-profile, and Figure 0 shows 
these approximations to the total luminosity as a func¬ 
tion of the Sersic index n. In the case of 2Rp when 
l/r/(Rp) = 0.2, one can see that profiles with n=10 will 
have their luminosities under-estimated by 44.7% and 
those with n=4 by only 17.1%. The situation is consid¬ 
erably worse when using 3Rp and l/r/(Rp) = 0.5. A 
prescription to correct for the missing light, beyond 
one’s chosen aperture, is detailed in Graham et al. 
(2005). 



Figure 7: Flux ratio, as a function of light-profile 
shape n, between the total luminosity L to t and the 
Petrosian luminosity Tp e t inside (i) twice the ra¬ 
dius i?p where l/?y(i?p) = 0.2 (solid curve) and (ii) 
thrice the radius Rp where l/ij(Rp) = 0.5 (dotted 
curve). 


This is a particular clever quantity because if ev¬ 
ery galaxy had the same stellar distribution, such as 
an l ? 1 / 4 profile, then a radius where the ^-function 
equalled some pre-defined, constant value would cor¬ 
respond to the same number of R e for every galaxy. 
Moreover, such measurements are unaffected by such 
things as exposure-depth, galactic dust, and cosmo¬ 
logical redshift dimming because they affect both sur¬ 
face brightness terms in Equation equally. Even 
though it is possible to measure the Petrosian radius 
without ever assuming or specifying an underlying light- 
profile model, the actual form of the stellar distribution 
is implicitly incorporated into the Petrosian function 
and so cannot be ignored (as Figure El reveals). 

It turns out the Petrosian function is equal to 


(29) 


2.6 Kron magnitudes 


Kron (1980) presented the following luminosity-weighted 
radius, R i, which defines the ‘first moment’ of an im¬ 
age 


Ri(R) 


2n /* I(x)x 2 dx 
27r f R I(x)xdx 


(31) 


He argued that an aperture of radius twice Ri, when 
Ri is obtained by integrating to a radius R that is 1% 
of the sky flux, contains more than ~90% of an object’s 
total light, making it a useful tool for estimating an 
object’s flux. 

It is worth noting that considerable confusion ex¬ 
ists in the literature in regard to the definition of Ri. 
To help avoid ambiguity, we point out that g(x ) in 
Kron’s (1980) original equation refers to xl(x), where 


V(R) = 2 /a(R) 
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x is the radius and I(x) the intensity profile. Infante 
(1987) followed this notation, but confusingly a typo 
appears immediately after his Equation 3 where he 
has written g(x) ~ I(x) instead of g(x) ~ xl(x). Fur¬ 
thermore, Equation 3 of Bertin & Arnouts (1996) is 
given as Ri = )P R/(R) / y) J(R) where the summa¬ 
tion is over a two-dimensional aperture rather than a 
one-dimensional light-profile. In the latter case, one 
would have Ri = J/ R 2 I(R)/ £ RJ(R). 

Using a Sersic intensity profile, and substituting in 
x — b(R/Re) 1 /" , the numerator can be written as 


27rn7 e e 6 7?e7(3n, x)/b 3n . 


Using Equation <n for the denominator, which is sim¬ 
ply the enclosed luminosity, Equation m simplifies 
to 


Ri(x, n) 


Re 7(3 n,x) 
b n 7(2n, x) 


(32) 


The use of ‘Kron radii’ to determine ‘Kron magni¬ 
tudes’ has proved very popular, and SExtractor (Bertin 
& Arnouts 1996) obtains its magnitudes using aper¬ 
tures that are 2.5 times Ri. Recently, however, it has 
been reported that such an approach may, in some in¬ 
stances, be missing up to half of a galaxy’s light (Bern¬ 
stein, Freedman, & Madore 2002; Benitez et al. 2004). 
If the total light is understood to be that coming from 
the integration to infinity of a light-profile, then what 
is important is not the sky-level or isophotal-level one 
reaches, but the number of effective radii that have 
been sampled. 

For a range of light-profiles shapes n, Figure © 
shows the value of Ri (in units of R e ) as a function of 
the number of effective radii to which Equation m 
has been integrated. Given that one usually only mea¬ 
sures a light-profile out to 3-4 R e at best, one can 
see that only for light-profiles with n less than about 
1 will one come close to the asymptotic value of Ri 
(i.e, the value obtained if the profile was integrated 
to infinity). Table © shows these asymptotic values 
of Ri as a function of n, and the magnitude enclosed 
within 2Ri and 2.5Ri. This is, however, largely aca¬ 
demic because observationally derived values of Ri will 
be smaller than those given in Table Q, at least for 
light-profiles with values of n greater than about 1. 


Table 1: Theoretical Kron Radii and Magnitudes 


Sersic n 

Ri 

[*e] 

L(< 2Ri) 

% 

L(< 2.5Ri) 

% 

0.5 

1.06 

95.7 

99.3 

1.0 

1.19 

90.8 

96.0 

2.0 

1.48 

87.5 

92.2 

3.0 

1.84 

86.9 

90.8 

4.0 

2.29 

87.0 

90.4 

5.0 

2.84 

87.5 

90.5 

6.0 

3.53 

88.1 

90.7 

7.0 

4.38 

88.7 

91.0 

8.0 

5.44 

89.3 

91.4 

9.0 

6.76 

90.0 

91.9 

10.0 

8.39 

90.6 

92.3 



Figure 8: Kron radii R\, as obtained from Equa¬ 
tion (El), are shown as a function of the radius R 
to which the integration was performed. Values of 
n range from 0.5, 1, 2, 3,... 10. 


From Figure © one can see, for example, that 
an f? 1 / 4 profile integrated to 4R e will have Ri=1.09R e 
rather than the asymptotic value of 2.29 R e . Now 2.5 x 
1.097?e encloses 76.6% of the object’s light rather than 
90.4% (see Table©. This is illustrated in Figure © 
where one can see when and how Kron magnitudes fail 
to represent the total light of an object. This short¬ 
coming is worse when dealing with shallow images and 
with highly concentrated systems having large values 
of n (brightest cluster galaxies are known to have Sersic 
indices around 10 or greater, Graham et al. 1996). 



Ri/R e 


Figure 9: Kron luminosity within 2.5-Ri, nor¬ 
malised against the total luminosity, as a function 
of how many effective radii R± corresponds to. Val¬ 
ues of n range from 0.5, 1, 2, 3,... 10. 


To provide a better idea of the flux fraction rep¬ 
resented by Kron magnitudes, and one which is more 
comparable with Figure ©, Figure m shows this 
fraction as a function of light-profile shape n. The 
different curves result from integrating Equation m 
to different numbers of effective radii in order to obtain 
RI. If n= 4, for example, but one only integrates out 
to lR e (where R e is again understood to be the true, 
intrinsic value rather than the observed value), then 
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the value of Ri is 0.41R e and the enclosed flux within 
2.5-Ri is only 50.7%. If an n=10 profile is integrated 
to only 1-Re, then Ri=0.30R e and the enclosed flux is 
only 45.0% within 2.5Ri. It is therefore easy to un¬ 
derstand why people have reported Kron magnitudes 
as having missed half of an object’s light. 


where Rb is the break-radius separating the inner power- 
law having logarithmic slope 7 from the outer Sersic 
function. The intensity R at the break-radius Rb can 
be evaluated from the expression 

/' = 7 6 2- (7/q) exp [fc(2 1 /Q R 6 /R e ) 1/n ] . (34) 



n 


Figure 10: Kron luminosity within 2.5Ri, nor¬ 
malised against the total luminosity, as a function 
of the underlying light-profile shape n. The dif¬ 
ferent curves arise from the different values of Rl 
obtained by integrating Equation (ITTT1) to (i) lR e , 
(ii) 2R e , (iii) 4R e , and (iv) infinity. 


2.7 The core—Sersic model 

Due to the presence of partially depleted cores in lumi¬ 
nous (Mb < —20.5 mag, Ho = 70 km s _1 Mpc^ 1 ) el¬ 
liptical galaxies 6 , a ‘core-Sersic’ model (Graham et al. 
2003a,b) has been developed in order to describe and 
connect the nuclear (typically less than a few hundred 
parsecs) and the remaining outer stellar distribution. 
This model consists of an inner power-law and an outer 
Sersic function, and has proven essential for model¬ 
ing the HST-resolved light-profiles of luminous early- 
type galaxies (Trujillo et al. 2004). As suggested in 
Graham et al. (2003b), the lower- luminosity ‘power- 
law’ galaxies have been shown to be described by the 
Sersic model over their entire radial extent (Trujillo et 
al. 2004). 

Although the reader is referred to the above pa¬ 
pers, especially the Appendix of Trujillo et al. (2004), 
the core-Sersic model is given as 7 : 




7/0 


exp{-6[(R Q + Rn/R:] 1/( “ n) } 
(33) 


6 The luminous ‘core galaxies’ likely correspond to the 
‘bright’ family of galaxies identified in Capaccioli, Caon, & 
D’Onofrio et al. (1992) and Caon et al. (1993). They tend to 
have boxy rather than disky isophotes (Nieto et al. 1991), 
and Sersic indices greater than ~4. They are commonly 
understood to be the product of (elliptical) galaxy mergers 
(e.g., Capaccioli, Caon, & D’Onofrio 1992, 1994; Graham 
2004) 

7 The a and 7 terms shown here should not be confused 
with those given earlier in the paper, they are different 
quantities. 


The final parameter, a, controls the sharpness of the 
transition between the inner (power-law) and outer 
(Sersic) regimes — higher values of a indicating sharper 
transitions. In practice (e.g., Fig-nre fTTIl we find that a 
sharp transition is adequate and recommend setting a 
to a suitably large constant value (typically anything 
greater than 3 is fine), leaving one with a 5-parameter 
core-Sersic model. 



radius [arcsec] 


Figure 11: Major-axis, R-band light-profile of 
NGC 3348. The solid line is the best fitting core- 
Sersic model while the dashed line is the best¬ 
fitting Sersic model to the data beyond the HST- 
resolved break radius (Graham et al. 2003a,b; Tru¬ 
jillo et al. 2004). 


2.8 Deprojected quantities and dy¬ 
namical terms 

Ciotti (1991) provides an exact, numerical expression 
for the deprojected light-profile of the R 1//n model, 
that is, the luminosity density profile. He additionally 
provides numerical expressions for the gravitational 
potential and also the spatial and line-of-sight veloc¬ 
ity dispersions. These must however be solved by inte¬ 
gration, which means they require considerably more 
computer time than analytical expressions. Ciotti does 
however provide analytical expressions for the behav¬ 
ior of the above expressions at both small and large 
radii. Luminosity-weighted aperture velocity disper¬ 
sions have been used in Ciotti, Lanzoni, & Renzini 
(1996), and also in Graham & Colless (1997) where 
the radial profiles are shown for different values of 
the Sersic index n. Ciotti (1991) additionally provides 
expressions for the distribution function and the nor¬ 
malised differential energy distribution. 
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An exact, analytical expression for the density pro¬ 
file was finally discovered a couple of years ago and is 
given in Mazure & Capelato (2002). It involves the 
use of somewhat complicated Meijer G functions. For 
those interested in a more simple, analytical approx¬ 
imation, an accurate expression is given in Prugniel 
& Simien (1997), which is developed slightly in Lima 
Nieto, Gerbal, & Marquez (1999) and Trujillo et al. 
( 2002 ). 

Mazure & Capelato (2002) also provide exact an¬ 
alytical expressions for the mass, gravitational poten¬ 
tial, total energy, and the central velocity dispersion. 
For modellers interested in fast-to-compute, analyti¬ 
cal approximations for not only the density profile but 
also the potential and force, such expressions, which 
additionally include optional power-law cores, can be 
found elsewhere (B. Terzic & A.W. Graham, in prepa¬ 
ration) . 

3 Sersic magnitudes 

In this article we have compiled and developed equa¬ 
tions pertaining to the Sersic profile in a purely an¬ 
alytical manner. To mention just one of many im¬ 
portant uses of the Sersic profile is its potential for 
deriving accurate total magnitudes. This need is mo¬ 
tivated by a growing community-wide awareness of the 
complex nature of galaxy photometry, and in partic¬ 
ular the large amounts of flux which can be missed 
by isophotal, aperture, Petrosian, or Kron magnitudes 
(e.g., Figures HI & llOtl . Cross et al. (2004) recently 
compared APM, SuperCosmos, SDSS and MGC pho¬ 
tometry for several thousand galaxies and concluded 
that the photometric errors are mainly dominated by 
the systematics of the methodology rather than the 
random errors. The Sersic magnitude provides a logi¬ 
cal standard and is derived by evaluating Equation^] 
at R = oo given p e , n, and R e derived from a fit to the 
available light-profile. 

In practice various “facts-of-life” issues remain; these 
are not specific problems to the Sersic model, but generic 
to studies of galaxy photometry. The most obvious 
ones are: the smoothing effect of the point-spread func¬ 
tion (PSF); profile truncation; multiple component sys¬ 
tems; dust and asymmetric profiles. All of these can 
act to exacerbate or ameliorate the amount of missing 
flux. While a detailed discussion of these issues is be¬ 
yond the scope of this paper we provide some suitable 
references below. 

The smearing effect of the PSF will cause the ob¬ 
served profile to tend to ti ~ 0.5 (i.e., Gaussian), this 
is dealt with straightforwardly by incorporating PSF 
convolution into the model fitting process, e.g., An- 
dredakis et al. (1995; their Equation 10) and Trujillo 
et al. (2001a,b). Disk truncation is trickier (see the 
reviews by van der Kruit (2001) and Pohlen et al. 
2004) and is assumed to be related to the minimum gas 
density required for star-formation (Kennicutt 1989; 
Kregel & van der Kruit 2004). Initially truncation 
was reported to occur at around 4 scale-lengths for ex¬ 
ponential disks (van der Kruit 1979; van der Kruit & 
Searle 1981). More recent studies have argued that the 


truncation is better described by a broken exponential 
fit (e.g., de Grijs et al. 2001; Pohlen et al. 2002). Oth¬ 
ers argue that truncation may actually be a manifes¬ 
tation of poor background modelling or simply due to 
intrinsic variations in the disk (Narayan & Jog 2003). 
Certainly some recent studies find no discernible trun¬ 
cation to extremely faint isophotal limits; for example, 
NGC 300 is a pure exponential out to 10 scale-lengths 
(Bland-Hawthorn et al. 2005). The net result is that in 
practice it is not clear exactly how far out one should 
integrate the Sersic profile for disk galaxies. As dis¬ 
cussed in Section ESI there is no evidence for trunca¬ 
tion in the elliptical galaxy population. From Figure^ 
we see that this issue is far more significant for high 
Sersic index systems. A new generation of publicly 
available 2D fitting codes, in particular GIM2D (Mar- 
leau & Simard 1998), GALFIT (Peng et al. 2002), and 
BUDDA (de Souza, Gadotti, & dos Anjos 2004), can 
routinely deal with multiple-component profiles. Dust 
opacity (Disney, Davies, & Phillipps 1989, Davies et 
al. 1993) can also lead to changes in the light-profile 
because of the more centrally concentrated dust distri¬ 
bution. Modelling opacity is non-trivial as there are 
strong inclination dependencies (Choloniewski 1991; 
Jones, Davies & Trewhella 1996; Graham 2001b) how¬ 
ever models are being developed to provide detailed 
corrections (e.g., Pierini et al. 2004; Tuffs et al. 2004). 
From the dust attenuation studies of Calzetti (2001, 
and references therein) and many others, dust issues 
can be minimised via structural analysis at near-IR 
wavelengths. Noil-biaxial asymmetry of galaxy images 
can be readily identified via the use of the ‘asymme¬ 
try’ measures (e.g., Schade et al. 1995; Rix & Zaritsky 
1995; Conselice 1997; Kornreich, Haynes, & Lovelace 
1998). 

Acknowledgments 

We kindly thank Valeria Coenda for faxing us a copy of 
Sersic’s 1963 article, referenced in his 1968 Atlas. We 
are additionally grateful to Steve Phillipps for referee¬ 
ing this work, and Massimo Capaccioli for his helpful 
and rapid response to our request for information and 
references. 

References 

[] Abraham, R.G., Valdes, F., Yee, H.K.C., & van den 
Bergh, S. 1994, ApJ, 432, 75 

Q Andredakis, Y.C., Peletier, R.F., & Balcells, M. 
1995, MNRAS, 275, 874 

[] Benitez, N., et al. 2004, ApJS, 150, 1 

[] Bernstein, R.A., Freedman, W.L., & Madore, B.F. 
2002, 571, 107 

(1 Bershady, M.A., Jangren, A., & Conselice, C.J. 
2000, AJ, 119, 2645 

[] Bertin, E., & Arnouts, S. 1996, A&AS, 117, 393 


10 

Publications of the Astronomical Society of Australia 

[] Bland-Hawthorn, J., Vlajic, M., Freeman, K.C., & 
Drain, B.T. 2005, AJ, submitted 

P de Vaucouleurs, G. 1977, in Evolution of Galaxies 
and Stellar populations, eds. R. Larson & B. Tinsley 

[] Blanton, M.R., et al. 2001, AJ, 121, 2358 

(New Haven: Yale University Observatory) p.43 

[] Blanton, M.R., et al. 2003, ApJ, 594, 186 

P Disney, M.J., Davies, J.I., & Phillipps, S. 1989, MN¬ 
RAS, 239, 939 

0 Calzetti, D. 2001, PASP, 113, 1449 

P Djorgovski, S., & Spinrad, H. 1981, ApJ, 251, 417 

[] Caon, N., Capaccioli, M., & D’Onofrio, M. 1993, 
MNRAS, 265, 1013 

P Djorgovski, S., Spinrad, H., & Marr, J. 1984, in 

New Aspects of Galaxy Photometry, J.L.Nieto (ed.), 

[] Caon, N., Capaccioli, M., D’Onofrio, M., & Longo, 

G. 1994, A&A, 286, 39 

Springer-Verlag 

P Doi, M., Fukugita, M., & Okamura, S. 1993, MN¬ 

0 Capaccioli, M. 1989, in The World of Galaxies, 
ed. H. G. Corwin, L. Bottinelli (Berlin: Springer- 
Verlag), 208 

RAS, 264, 832 

P D’Onofrio, M. 2001, MNRAS, 326, 1517 

[] Capaccioli, M., Caon, N., & D’Onofrio, M. 1992, 
MNRAS, 259, 323 

P D’Onofrio, M., Capaccioli, M., & Caon, N. 1994, 

MNRAS, 271, 523 

P Capaccioli, M., Caon, N., & D’Onofrio, M. 1994, 
Mem. della Societa Astronomia Italiana, 65, 919 

P Driver, S.P., Phillipps, S., Davies, J.L, Morgan, I., 

& Phillipps, S. 1994, MNRAS, 266, 155 

P Capaccioli, M., & de Vaucouleurs, G. 1983, ApJS, 

52, 465 

P Ellis, G.F.R., & Perry, J.J. 1979, MNRAS, 187, 357 

P Ferrari, F., Dottori, H., Caon, N., Nobrega, A., & 

P Cellone, S.A., Forte, J.C., & Geisler, D. 1994, ApJS, 

93, 397 

Pavani, D.B. 2004, MNRAS, 347, 824 

P Fraser, C.W. 1972, The Observatory, 92, 51 

P Choloniewski, J. 1991, MNRAS, 250, 486 

P Graham, A.W. 2001a, AJ, 121, 820 

P Ciotti, L. 1991, A&A, 249, 99 

P Graham, A.W. 2001b, MNRAS, 326, 543 

P Ciotti, L., & Bertin, G. 1999, 352, 447 

P Graham, A.W. 2004, ApJ, 613, L33 

P Ciotti, L., Lanzoni, B., & Renzini, A. 1996, MN¬ 
RAS, 282, 1 

P Graham, A.W., & Colless, M.M. 1997, MNRAS, 

287, 221 

P Conselice, C.J. PASP, 1997, 109, 1251 

P Graham, A.W., Erwin, P., Trujillo, I., & 

P Cross, N.J.G., Driver, S.P., Liske, J., Lemon, D.J., 
Peacock, J.A., Cole, S., Norberg, P., & Sutherland, 
W.J. 2004, MNRAS, 349, 576 

Asensio-Ramos, A. 2003a, in Carnegie Ob¬ 
servatories Astrophysics Series, Vol. 1: Co¬ 

evolution of Black Holes and Galaxies (ed. 

L. C. Ho (Pasadena: Carnegie Observatories, 

P Dalcanton, J.J., & Bernstein, R.A. 2000, AJ, 120, 
203 

http: //www.ociw.edu / ociw/symposia / series / symposiuml /proceeding 

P Graham, A.W., Erwin, P., Trujillo, I., & Asensio- 

[I Davies, J.I., Phillipps, S., Boyce, P.J., & Disney, 
M.J. 1993, MNRAS, 260, 491 

Ramos, A. 2003b, AJ, 125, 2951 

P Graham, A.W., & Guzman, R. 2003, AJ, 125, 2936 

P Davies, J.I., Phillipps, S., Cawson, M.G.M., Disney, 
M.J., & Kibblewhite, E.J. 1988, MNRAS, 232, 239 

P Graham, A.W., Lauer, T., Colless, M.M., & Post¬ 
man, M. 1996, ApJ, 465, 534 

P de Grijs, R., Kregel, M., & Wesson, H. 2001, MN¬ 
RAS, 324, 1074 

P Graham, A.W., Trujillo, I., & Caon, N., 2001, AJ, 

122, 1707 

P de Souza, R.E., Gadotti, D.A., & dos Anjos, S. 2004, 
ApJS, 153, 411 

P Graham, A.W., et al. 2005, AJ, submitted 

P de Vaucouleurs, G. 1948, Ann. Astrophys., 11, 247 

P Gunn, J.E., & Oke, J.B. 1975, ApJ, 195, 255 

P de Vaucouleurs, G. 1956, Mem. Commonwealth 
Obs., 3, No. 13 

P Hoessel, J.G. 1980, ApJ, 241, 493 

P Infante, I. 1987, A&A, 183, 177 

P de Vaucouleurs, G. 1959, Handbuch der Physik, 53, 
275 & 311 

P Iodice, E., D’Onofrio, M., Capaccioli, M. 1997, ASP 

Conf. Ser., 116, 841 


www.publish.csiro.au/journals/pasa 


11 


[] Iodice, E., D’Onofrio, M., Capaccioli, M. 1999, ASP 
Conf. Ser., 176, 402 

[] Jones, H., Davies, J.I., & Trewhella, M., 1996, MN- 
RAS, 283, 316 

[] Karachentseva, V.E., Prugniel, P., Vennik, J., 
Richter, G.M., Thuan, T.X., & Martin, J.M. 1996, 
A&ASS, 117, 343 

[] Kennicutt, R.C. 1989, ApJ, 344, 685 

[] Khosroshahi, H.G., Wadadekar, Y., & Kembhavi, 

A. 2000, ApJ, 533, 162 

[] Kjaergaard, P., Jorgensen, I, & Moles, M. 1993, ApJ, 
418, 617 

[] Kornreich, D.A., Haynes, M.P., & Lovelace, R.V.E. 
1998, AJ, 116, 2154' 

[] Kregel, M., & van der Kruit, P.C. 2004, MNRAS, 
355, 143 

[] Kron, R. 1980, ApJS, 43, 305 

[] Lima Nieto, G.B., Gerbal, D., & Marquez, I. 1999, 
MNRAS, 309, 481 

[] MacArthur, L.A., Courteau, S., & Holtzman, J.A. 
2003, ApJ, 582, 689 

[] Marleau, F.R., & Simard, L. 1998, ApJ, 507, 585 

[] Mazure, A., & Capelato, H.V. 2002, A&A, 383, 384 

[] Mollenhoff, C., & Heidt, J. A&A, 368, 16 

[] Morgan, W.W. 1958, PASP, 70, 364 

[] Morgan, W.W. 1959, PASP, 71, 394 

[] Morgan, W.W. 1962, ApJ, 135, 1 

[] Narayan, C.A., & Jog, C.J. 2003, A&A, 407, L59 

[] Nieto, J.-L., Bender, R., & Surma, P. 1991, A&A, 
244, L37 

[] Okamura, S., Kodaira, K., & Watanabe, M. 1984, 
ApJ, 280, 7 

[] Peng, C.Y., Ho, L.C., Impey, C.D., & Hans-Walter, 
W. 2002, AJ, 124, 266 

[] Petrosian, V. 1976, ApJ, 209, LI 

[] Pierini, D., Gordon, K.D., Witt, A.N., & Madsen, 
G.J. 2004, ApJ, 617, 1022 

[] Pohlen, M., Dettmar, R., Lutticke, R., & Aronica, 
G. 2002, A&A, 392, 807 

[] Pohlen M., Beckman, .J.E., Hiittemeister, S., 
Knapen, J.H., Erwin, P., & Dettmar, R.-J. 2004, in 
Penetrating Bars Through Masks of Cosmic Dust, 
eds., D.L. Block, I. Puerari, K.C. Freeman, R. 
Groess, E.K. Block, A&SS, 319, 713 

[] Prugniel, P., & Simien, F. 1997, A&A, 321, 111 


[] Rix, H.-W., & Zaritsky, D. 1995, ApJ, 447, 82 

[] Rest, A., et al. 2001, AJ, 121, 2431 

[] Sandage, A., & Perelmuter, J.-M. 1990, ApJ, 350, 
481 

[] Schade, D., Lilly, S.J., Crampton, D., Hammer, F., 
Le Fevre, O., & Tresse, L. 1995, ApJ, 451, LI 

[] Seigar, M.S., & James, P.A. 1998, MNRAS, 299, 672 

[] Sersic, J.-L. 1963, Boletin de la Asociacion Ar¬ 
gentina de Astronomia, vol.6, p.41 

[] Sersic, J.-L. 1968, Atlas de Galaxias Australes (Cor¬ 
doba: Observatorio Astronomico) 

0 Shen, S., et al. 2003, MNRAS, 343, 978 

[] Trujillo, I., Aguerri, .J.A.L., Cepa, J., & Gutierrez, 
C.M. 2001a, MNRAS, 321, 269 

[] Trujillo, I., Aguerri, .J.A.L., Cepa, J., & Gutierrez, 
C.M. 2001b, MNRAS, 328, 977 

[] Trujillo, I., Asensio Ramos, A., Rubino-Martin, 
J.A., Graham, A.W., Aguerri, J.A.L., Cepa, J., & 
Gutierrez, C.M. 2002, MNRAS, 333, 510 

[] Trujillo, I., Erwin, P., Asensio Ramos, A., & Gra¬ 
ham, A.W. 2004, AJ, 127, 1917 

[] Trujillo, I., Graham, A.W., & Caon, N. 2001c, MN¬ 
RAS, 326, 869 

0 Tuffs, R., Popescu, C.C., Volk, H.J., Kylafis, N.D., 
& Dopita, M.A. 2004, A&A, 419, 821 

[] van der Kruit, P.C. 1979, A&A, 38, 15 

[] van der Kruit, P.C. 2001, in Galaxy Disks and Disk 
Galaxies, G. Jose, S.J. Funes, E.M. Corsini, ASP 
Conf. Ser., 230, 119 

[] van der Kruit, P.C., & Searle, L. 1981, A&A, 95, 
105 

[] Vennik, J., Hopp, U., Kovachev, B., Kuhn, B., & 
Elsasser, H. 1996, A&ASS, 117, 261 

[] Vennik, J., & Richter, G.M. Astron. Nachr., 315, 
H3, 245 

[] Young, C.K., & Currie, M.J. 1994, MNRAS, 268, 
Lll 

[] Young, C.K., & Currie, M.J. 1995, MNRAS, 273, 
1141 



